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Abstract 

In the present paper we analyze the linear stability of a hierarchical 
size-structured population model where the vital rates (mortality, fertility 
and growth rate) depend both on size and a general functional of the 
population density ("environment"). We derive regularity properties of 
the governing linear semigroup, implying that linear stability is governed 
by a dominant real eigenvalue of the semigroup generator, which arises as 
a zero of an associated characteristic function. In the special case where 
neither the growth rate nor the mortality depend on the environment, 
we explicitly calculate the characteristic function and use it to formulate 
simple conditions for the linear stability of population equilibria. In the 
general case we derive a dissipativity condition for the linear semigroup, 
thereby characterizing exponential stability of the steady state. 

Keywords: Hierarchical size-structured populations; Semigroup meth- 
ods; Spectral analysis; Principle of linear stability 



1. Introduction 

In the last three decades nonlinear age- and size-structured population mod- 
els have attracted a lot of interest both among theoretical biologists and applied 
mathematicians. Traditionally, structured population models have been formu- 
lated as partial differential equations for population densities. Starting with 
the seminal work |16j . researchers have been developing and analyzing various 
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physiologically structured population models. We refer here to the well-known 
monographs [TJ H2 [25] . 

Diekmann et al. have been developing a general mathematical framework 
for modeling structured populations, see for example [9j [10] . One of their most 
important recent results is that the qualitative behavior of nonlinear physio- 
logically structured population models can be studied by means of linearization 
[3 [8] . In other words, they have proven for a very general class of physiologically 
structured population models that the nonlinear stability/instability of a pop- 
ulation equilibrium is completely determined by its linear stability /instability. 
Such a fundamental result, often referred to as "the Principle of Linearized Sta- 
bility" , has been shown previously for several concrete age- and size-structured 
models [M [H [H [H \M . 

Following the lead of [23] and [25], we successfully applied linear semi- 
group methods to formulate biologically interpretable conditions for the lin- 
ear stability/instability of equilibria of several structured population models 
[T3"l ITU [T5] . In these problems the vital rates depend on size or age and on 
the total population size, in general. Hence it is assumed that any effect of 
intraspecific competition on individual behavior is primarily due to a change 
in population size and every individual in the population can influence the vi- 
tal rates of other individuals, a scenario commonly referred to as "scramble 
competition" . 

In other scenarios competition among individuals is based upon some hier- 
archy in the population which is often related to the size of individuals. In this 
case the nonlinearity (environmental feedback) in the model is incorporated 
through infinite dimensional interaction variables. A simple example for this 
situation is given by a forest consisting of tree individuals in which the height of 
a tree determines its rank in the population [5T] . Taller individuals have higher 
efficiency when competing for resources such as light, while individuals of lower 
rank cannot affect the vital rates of individuals of higher rank. This scenario 
is the so called "contest competition". Both discrete time and continuous hi- 
erarchical structured models have been developed, see [6] and the references 
therein. 

Of interest in this work is the stability analysis of population equilibria by 
means of linearization of a continuous quasilincar size-structured model, recently 
discussed in [2j. In this model the density evolution of individuals of size s is 
assumed to be governed by the following quasilinear partial differential equation 



defined for < s < m < oo and t > 0. The density of zero (or minimal) size 
individuals is given by the nonlocal boundary condition 



Jo 

The quantity m denotes the maximum size of individuals. The initial condition 
takes the form 



u t (s, t) + (j(s, Q{s, t)) u(s, t)) s + fi(s, Q(s, t)) u(s, t) = 0, 



(1.1) 




(1.2) 



u(s, 0) = uq(s), s€[0,m}. 



(1.3) 
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Here /?, /i and 7 denote the fertility, mortality and growth rate of individuals, 
respectively. We assume that these vital rates depend on the individual size s 
and on the environment 

Q(s,t) — a / w(r]) u(rj,t) drj + / w(rf)u(tj,t)drj, < s < m, t > 0. 

(1.4) 

The constant a is a parameter in [0, 1] measuring the degree of hierarchy in the 
population, while the function w represents a positive weight. For example in 
case of a tree population where taller individuals overshadow smaller individuals 
|21j the vital rates of an individual of size s are reasonably assumed to depend 
on the cumulative leaf area of individuals of size s or larger, modeled by the 
function 

w(r))u(r),t)dr). (1.5) 

Here w is an appropriately chosen weight function. Hence in this situation the 
parameter a would be 0. The case a = 1 (which represents scramble competi- 
tion) has been treated in detail in [IB] . 

We impose the following regularity conditions on the model ingredients: 

• fi = n(s,Q) e C([0, m];^ [Q,c»)), n > 

• 7 = 70, Q) e C x {[Q,m]\C x % 00)) n C([0,m]; C 2 [0, 00)), 7 > 

• = P(s,Q) eC([0,m];C 1 [0,oo)),/3>0 

• w = w(s) e (^([O, m]), w > 0. 

These assumptions are tailored toward the linear analysis of this work. They 
might, however, not suffice to guarantee the existence and uniqueness of solu- 
tions of Eqs. (|1.1|) - (|1.4|) . Well-posedness of structured partial differential equa- 
tion models with infinite dimensional environmental feedback variables is in 
general an open question, although conditions for the global existence of weak 
solutions in the case discussed here are given in [2j . It has recently been shown 
[2 HT] that the population model (|l.ip - (|1.4p may exhibit a more complicated 
dynamical behavior than the simple size-structured model of scramble competi- 
tion. In particular, in [T] it was demonstrated both analytically and numerically 
that a singular solution of (|1.1|) - (|1.4|) containing a Dirac delta mass component 
can emerge if the growth rate 7 is not a decreasing function of the environment 

Q- 

For a more realistic description of real populations in a specific setting, one 
will have to modify the assumptions on the vital rates above. For example, one 
would possibly demand that lim /z(s, .) = 00, thus modeling a gradual rather 

s — >?n 

than instantaneous reduction in the numbers of individuals reaching maximum 
size to. 
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The size-structured model (ll.ll) - (|1.4p is often considered (see [TJ [2] ) with a 
boundary condition of the form 



7(0,Q(0,i))it(0,t) = C(t)+ / /3(s,Q(s,t))u{s,t)ds, t > 0. (1.6) 



In (|1.2[) we have taken C = and incorporated the growth rate 7(0, Q(0, t)) on 
the left of (|1.6p in the birth rate j3(s, Q(s, t)) on the right of (|1.2[) . assuming that 
zero size individuals grow instantaneously. This assumption seems reasonable, 
for example in case of a forest population. It is then clear that the two boundary 
conditions are equivalent (in the case C = treated here). We have observed, 
however, that (|1.2p is better suited for analytical work [13l HH [15] . As recent 
results indicate P21 [15], the introduction of a positive inflow C may have a 
significant influence on the linearized dynamical behavior of (|1.1|) - (|1.4|) . A 
comprehensive study of the effects of a positive inflow in hierarchical populations 
is left for future work. 

The study of hierarchical models in the literature is largely based on a de- 
coupling of the total population quantity from the governing equations and a 
transformation of the nonlocal partial differential equation into a local one 
0H1 [21]. This technique allows to prove well-posedness and to study the asymp- 
totic behavior of solutions by means of ODE methods. For Eqs. (|1.1|| - (|1.4|) this 
transformation fails since the vital rates depend on both size s and on the en- 
vironment Q(s,t). Therefore it seems unavoidable to study the original partial 
differential equation (jl.ip with the nonlocal integral boundary condition (|1.2[1 
directly. This approach is based on a linearization of the governing equations 
about steady state [H EH HI E2 H3 123 . While Sections 3 through 5 exploit 
spectral theoretic and structural properties of the governing linear semigroup 
extending related results in [13] [14] [15] , Section 6 gives a new characterization 
of asymptotic stability of the semigroup in terms of a dissipativity criterion. 
This idea was previously introduced and employed in |17[ 118] for elongational 
flow problems. 



Eqs. (|l.ip - (|1.4p have obviously the trivial solution it* = 0. Realistically we 
also expect additional positive (continuously differentiable) solutions it* > 0. In 
the following we formulate a necessary condition for the existence of a positive 
equilibrium solution of problem l|l.l|) - (|1.4|) . 

Proposition 2.1 Ifu* is a positive stationary solution of problem (|l.l|) - (|1.4jl . 

then the function Q*, defined by 



■m 



2. The linearized system 




(2.1) 



satisfies the equation 



R{Q*) = 1, 



(2.2) 
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where R : C([0, m]) — > R is t/ie inherent net reproduction rate 

R(Q)= P{s,Q(s))n(s,Q)ds (2.3) 
Jo 

and </ie operator it is given for < s < m and Q G C([0, m]) &?/ 

0) ^ exp (- f 4^ drl . (2.4) 



l(s,Q(s)) r I 7o 7(r,Q(r)) 

Proof. For a positive stationary solution let be given by (|2.1[) . Since 
any stationary solution satisfies 

u*(s) = u*(0) 7r(s, (2.5) 

we obtain Eq. (|2.2p when imposing the boundary condition (jl.2p . 



Given any stationary solution it* in C ([0, m]), we linearize the governing 
equations by introducing the infinitesimal perturbation v = v(s,t) and making 
the ansatz u — v + u* . After inserting this expression in the governing equations 
and omitting all nonlinear terms, we obtain the linearized problem 

V t (s,t) + 7* (s) v s (s, t) + p* (s) v(s, t) + a„ (s) V(s, t) = 0, (2.6) 

r rn r rn 

v(0,t)= (3(s,Q*(s))v(s,t)ds+ (3 Q (s,Q*(s))u4s)V(s,t)ds, (2.7) 



where we have set 

V(s,t)=a / w(rf) v(r),t) drj + / w{rf)v{r],t)dr], (2-8) 

7*(*)=7(«,Q.(*)). ( 2 -9) 
p,(s) =|i(s,Q*(a))+7,(«,Q.(«))+2(a-l)«;(s)7g(s > Q,(a))«,(a), (2.10) 

= Mq( s >Q*( s )) u *( s ) +7sq(s,Q*( s )) u *( s ) + 1q( s ,Q*( s )) u '*( s ) 

+ (a-l) w(s) 7qq(s, Q»(s)) u*(s) 2 . (2.11) 

We denote the Lebesgue space L (0,m) with its usual norm || • || by X and 
introduce the bounded linear functional A on X by 



A(v) = / /3(s,Q*(s))v(s)ds (2.12) 
Jo 

+ / Pq{s,Q*{s))u*(s) [a I w{rj) v(rj) dr\ + / w{rj) v(rj) dr\ I ds. 



Next we define the operators 

Av = -j(-,Q*)v s , Dom(A) = {v€W 1 < 1 (0,m)\v(0) = A(v)}, (2.13) 
Bv = -p*v on X, (2.14) 



Cv = — o~* [a / w(rj) v(rj) drj + / w(ji)v(ri) dr/ I on ^f. (2-15) 
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Then the linearized system (|2.6[) - (|2.7p can be cast in the form of an initial value 
problem for an ordinary differential equation on X 

j f v= (A + B + C) v, (2.16) 
together with the initial condition 

«(0) = wo. (2.17) 

In analogy to previously discussed size-structured population models [T51 [T^l [T5] . 
we can invoke the Desch-Schappacher Perturbation Theorem [11] to obtain the 
following result. 

Proposition 2.2 The operator A + B + C generates a strongly continuous semi- 
group {T(t)} t >o of bounded linear operators on X. 

The proof is a minor modification of parallel results given in [T31 HH [H] and 
has therefore been omitted. 



3. Spectral analysis and semigroup regularity 

Proposition 3.1 The spectrum of A + B + C can contain only isolated eigen- 
values of finite multiplicity. 

Proof. We prove that the resolvent operator of A + B + C is compact. Since 
B + C is a bounded perturbation of A, it suffices to show that the resolvent 
operator of A is compact. To this end, given / <E X, we find a unique solution 
v G Dom(.A) of the equation 

Xv-Av = f (3.1) 

in the form 

»« = ^ rw ( A( »» + I eArw TO))*) (3 - 2) 

if A € M. is sufficiently large. Here we define 

r« = f 4t dr >- ^ 

Jo 7*W 

Consequently, for A > large enough, the resolvent operator (XT — A)^ 1 exists 
and is bounded, mapping X = L 1 (0,to) into W /1,1 (0,to). Since W 1 ' 1 (Q, m) is 
compactly embedded in -^(0, m), the claim follows. 

□ 
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Theorem 3.2 The semigroup {T(i)} t >o, generated by the operator A + B + C , 
is eventually compact. Conseguently, the Spectral Mapping Theorem holds true, 
i.e. 

a(T(t)) = {0}Uc^p{a{A + B + C)t}, t > 0. (3.4) 

Moreover, the semigroup is spectrally determined, i.e. the growth rate to(T) of 
the semigroup and the spectral bound s(A + B + C) of its generator coincide. 

Proof. Since the operator C is compact, it is enough to prove the claim for 
the operator A + B. The differential equation 

^v = (A + B)v (3.5) 
dt 

corresponds to the partial differential equation 

v t (s,t) +7*(s)u s (s,t) + p*(s)v(s,t) = (3.6) 
together with the boundary condition (|2.7p . For to > let us introduce 

w(s)=«(a, *(«)), (3.7) 

where 

t(s) = t + T(s). (3.8) 

Then to satisfies the equation 

LO i {s) + £M uj{s) = Q, (3.9) 

7*( s ) 

hence 

w ( s ) = *o)) tt(s, Q*) exp I (1 - a) J - d-qj . 

(3.10) 

Thus for t — T(s) > we have 

v{s,t) = A(v{;t-T(s)))n{s,Q*) exp I (1 - a) J - ^ dr? I . 

(3.H) 

Therefore, noting the definition of A in (|2.12j) . we conclude that v is continuous 
in s and t if i > T(m) = maxo< s < m T(s). Consequently, Eq. (|3 . 1 1 1) in combi- 
nation with Eq. (|3.6[) implies that v is continuously differentiable if t > 2T(m). 
Hence the semigroup generated by A+B is differentiable for t > 2 T(m). Finally, 
since W 1,:L (0, to) is compactly embedded in L 1 (0, to), the semigroup is compact 
for t > 2T(m). The validity of the Spectral Mapping Theorem and the claim 
about the spectral determinacy of the semigroup follow, see [TT] . 

□ 

We conclude this section by formulating conditions for the positivity of the 
semigroup {T(t)} t >o- 
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Theorem 3.3 Suppose that 
er* < and (3.12) 

P(-,Q*)+w ^ /3 Q (ri,Q*(r)))u*(r))dri + a J /3q(v,Q*(v)) u *(v) d vj > 0- 

(3.13) 

TTien ifte semigroup {T(t)} t >o, generated by the operator A + B + C , is positive. 

Remark 3.4 Conditions (|3.12ll . (|3. 13|) are immediate generalizations of the 
corresponding positivity conditions given by Priifi in [23] for an age-structured 
scramble competition model. In general, if (3q = 0, condition (|3.13[) is trivially 
satisfied. Also, if the growth rate is independent of the environment (i.e. 7 = 
j(s)), condition (|3.12|) reduces to 

<0, sG[0,oo). (3.14) 

Hence in this case mortality is required to be a non-increasing function of the 
environment as well. 

Proof of Theorem 13.31 Since C is a positive operator by condition (I3.12|) . it 

suffices to prove the claim for the semigroup generated by A + B. Hence we 
assume that v satisfies Eq. (|3.6j) such that the boundary condition (|2.7|) and the 
initial condition v — vq 6 Dom(^4) for t = hold true. Let the function e* be 
given by 

e*(s)=exp 1(1 -a) / - -—^ drj \ . (3.15) 

V Jo 1*(V) J 



Then the function cj>, defined by 



^ ' ^ 7r(s, Q») e*(s) ' 3 6 



solves the problem 



M*.*)+7(*,Q*(*))0»(a.*)=O, (3.17) 
0(O,*)=A(^(.,t)7r(-,Q,)e,), (3.18) 

^.°)= f m = *»(')• ( 3 - 19 ) 

7r(s, <y*) e*(s) 

This boundary-initial value problem corresponds to the abstract initial value 
problem 

^ = Am <t>, m = 0o (3.20) 
at 

with the modified semigroup generator Am, defined by 

Ajrfcj) = — 7(-, Q*) 4> s on the domain 

Bom(A M ) = {</> G W 1 ' 1 (Q,m) | 0(0) = A (0tt(-, Q.) e*)} . (3.21) 
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For A > and g 6 £ (0, m), the resolvent equation 

\4>-A M <t> = 9 (3.22) 

has the implicit solution 

0(.s) = A (<M, Q.) e.) + f C^" 1 ™) ^ dr. (3.23) 

Applying A, we deduce the equation 

A /y e A(r W -r(.)) 9%) 

\Jo l{r,Q*[r)) 

l-A(e- xr ir(-,Q*)e*) 



if A is large enough. Condition (|3 . 1 3[) guarantees that A is a positive linear 
functional. Hence the solution 4>, given by Eq. (|3.23p . is nonnegative if g is 
nonncgative and A is sufficiently large, ft follows that the resolvent operator of 
Am (and consequently of A+B) is positive if A is large enough. This observation 
proves the claim. 

□ 

The positivity of the semigroup has far-reaching consequences. In particular, 
we obtain the following result from the theory of positive semigroups [llj . 

Corollary 3.5 Suppose that the semigroup {T(t)} t >o, generated by the operator 
A + B + C, is positive. Then the spectral bound s(A + B + C) € [-co, oo) satisfies 

s(A + B + C) = max {A e M | A is eigenvalue of A + B + C}. (3.25) 

Moreover, the spectrum of A + B + C is nonempty if and only if the spectral 
bound is finite. 



4. The characteristic equation 

In light of Theorem 13.21 the growth of the governing semigroup is deter- 
mined by the eigenvalues of its generator. Hence it is essential to determine the 
eigenvalues of A + B + C. The eigenvalue equation 

Aw - {A + B + C)v = (4.1) 

for A 6 C and nontrivial v is equivalent to the system 

v'(s) 1 4s)+v(s)(\ + p4s)) + V(s)a4s) = 0, (4.2) 

r 771 rTTl 

v(0)= p(s,Q*(s))v(s)ds+ f3 Q (s,Q*{s))u*( s )V(s)ds, (4.3) 
Jo Jo 
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where 

V(s) = a / w{rj) v{rj) d/q -\- I w(rj) v{rj) dr\ 

= (a — 1) / 10(77) 17(77) c?7/ + / w(r]) v(rj) dr\. (4-4) 



For the remainder of this section let us assume that a £ [0, 1). From (|4.4|) we 
obtain 

y'( s ) = ( a -l)io(s)u(s) and V"(s) = (a-1) (w'(s) v(s)+w(s) v'(s)). (4.5) 

Using the relations (|4.5p we can rewrite system (|4.2[) - (|4.3|) in terms of V and 
its derivatives as follows 

V"(s) + V'(s) + * - ^M) + V(s) (a - 1) «,(.) ^4 = 0. (4.6) 

V 7*( s ) w{s) J 7*(s) 

Eq. ()4.6j) is accompanied by boundary conditions of the form 

aV(0) = V(m), (4.7) 

V'(0)=w(0) r P( S ' Q *W V'(s)ds 
Jo w{s) 



+ w(Q)(a-l) / p Q (sMs))u*(s)V(s)ds. (4.8) 
Jo 

For A £ C, any solution V\(s) of the second order homogeneous ordinary differ- 
ential equation (|4.6j) can be written as 

V x (s)=c 1 V 1 (s,X) + c 2 V 2 (s,X), (4.9) 

where Vi(s,A) and V2(s,A) are any fixed, linearly independent solutions of 
Eq. (|4.6p and Ci, C2 are arbitrary constants. When imposing the boundary 
conditions (|4.7|) - (|4.8p . we obtain the conditions 

c 1 V{(0 > X) + c 2 Vl(0,\)=c 1 f ^f3{ S) Q*{s))Vl( S) \)ds (4.10) 

Jo w{s) 

+ c 2 r^p(s ) Q*(s))Vl(s ) \)ds (4.11) 
Jo w ( s ) 



+ ci / (a-l)u;(0)/3o(* s Q,(s))tt*(8)yi(fl,A)dfl (4.12) 
Jo 

r-ra 

+ c 2 / (a-lM0)/3 Q (s,Q»(s))w*(s)U 2 (s,A)ds (4.13) 



or in short 

ciJJi(A)+c 2 ff 2 (A) =0, (4.14) 

and 

ci a Vi(0, A) + c 2 a V 2 (0, A) = ci Vi(m, A) + c 2 V 2 (m, A), (4.15) 
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in short 

ci Ji(A)+c 2 J 2 (A) = 0. (4.16) 

Here the functions Hi, H2, Ji, and J2 represent the terms multiplying ci, 
C2, respectively. The homogeneous system (I4.14[) . (|4.16l) admits a nontrivial 
solution for c\ , ci if and only if A satisfies the equation 

#i(A) J 2 (A)-# 2 (A)Ji(A) =0. (4.17) 

This equation is the characteristic equation of the linearized system (|2.6[) - (|2.7p . 
Its zeros are the eigenvalues of the operator A+B+C, which completely describe 
the spectrum of A + B + C. 

The explicit information contained in the characteristic equation is, however, 
rather limited since linearly independent solutions of the second order differen- 
tial equation (|4.6p are in general not directly available, unless one resorts to 
numerical techniques. As we will see in the forthcoming section this problem 
can, however, be overcome in special cases of the model ingredients. 



5. A special case 

In this section we treat the special case when the mortality and growth rate 
are independent of the environment Q, i.e. 7q = = hq. Hence c* = and e* = 
1. In this situation we are able to determine the characteristic equation (|4. 1 7|) 
explicitly and to formulate simple conditions for the linear stability /instability 
of positive stationary solutions. In contrast to the preceding section we allow 
a G [0,1]. 

Theorem 5.1 Suppose a* = 0. Then a positive stationary solution is lin- 
early asymptotically stable if 

Pq(;Q*)<0, /3 q (-,Q*)#0 (5.1) 
and the positivity condition (|3.13jl holds true. If, however, 

Pq(-,Q*)>0, /3q(-,Q*)^0, (5.2) 
then it* is linearly unstable. 

Note that the instability part of the theorem does not require the positivity 
condition. 

Proof. We assume first that < a < 1. Then the general solution of (|4.6[) is 
found as 

V(s) = V(0) + V'(0) J* H(A, r) dr, (5.3) 

where we have set 

m . s def 7*(Q) / / S A + ^*(r) 1 

n (A, r) = —— exp <^ - / — — dr \ . (5.4) 

7*00 I Jo 7*W J 
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Imposing the boundary condition (|4.7[) on the solution (|5.3p , we obtain 

= 7(0) (l-a) + V'(0) j™ ^ n(A, s) ds, (5.5) 
while the boundary condition (|4.8p gives 

= 7(0) ft»(0) (1 - a) /3 Q ( S , Q,(«)) ds 

+ 7(0) U-jJ /5(s,Q*(s))n(A,s)ds (,-».(») 

+(l-a) / /3q(s,Q*(s))u*(s) / w(r) U(X, r) dr ds 



The linear system (|5.5[) - (l5.6p has a nontrivial solution (7(0), 7'(0)) if and only 
if A satisfies 



1 = / /3(s,Q*(s))H"(A,s)ds 
Jo 

+ (o-l) / f3 Q (s,Q*(s))u*{s) I w(r)n{X,r)drds (5.7) 



+ / f3 Q (s,Q*(s))u*(s)ds w(s)Tl(\,s)ds d = K{\). 
Jo Jo 

This equation corresponds to the characteristic equation (|4. 1T[) . If, however, 
a = 1, 7, defined by l|4.4j) . is constant. When we solve the problem ()4.2[) - (|4.3|) 
directly, we obtain again the condition 

K(\) = 1, (5.8) 



where K is given by (|5.7|) with a = 1. Hence (j5.7[) is the characteristic equa- 
tion for all < a. < 1. For the stability part, our assumptions guarantee that 
the positivity conditions (|3.12|) . (|3.13j) hold true. Therefore, to prove asymp- 
totic stability, it suffices to show that the characteristic equation (|5.7p has no 
nonnegative (real) solutions. To this end, we observe that 



K(0) =R(Q,) + / (3q(s, Q*(s)) u,(s) 
Jo 

X a / w(r) 7r(r, Q*) dr + / w(r) 7r (r, Q») dr I ds < 1 (5.9) 



by condition f)5 . 1 [) . Moreover, the positivity condition (|3 . 13[) yields that 



+ / /3 Q faQ*(r)KMdr (5.10) 

Jo 

+aw(s) / Pa{r,Q*{r))u*{r) dr ds < 0. 
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Consequently, K is monotone decreasing for A > 0. Hence the stability part is 
proven. The instability part of the theorem follows from the Intermediate Value 
Theorem since K(0) > 1 by ([5~2]) and lim K(\) = 0. 

A^oo 

□ 

Example 5.2 Let us consider an example where Theorem 15.11 yields asymp- 
totic stability. We choose 

m = 1, a — — , w = 1 (5-11) 

and let 

7 ( S ) = 1-I S , M =l, (5.12) 

480 ^ 
/9(*,Q)= 997(! + s )(3-2Q) if Q < ^ < s < 1, (5.13) 

where we assume that (5 extends to a continuously differentiable, non-negative 
function on [0, 1] x [0, oo). Then the corresponding problem has the stationary 
solution n,(s) = l-is with 

Q*(s) = j-^ + l< \ for0<*<l. (5.14) 
It is readily seen that 

/3q(s,Q*(s)) = -^(1 + s)<0, 0<s<1 (5.15) 

and that the positivity condition (|3.13[) reduces to 

s 3 s 2 3 s 5 , 
-^ + -r + -r + 7Tj>0, 0<s<l. 5.16 
24 4 4 24 

Since this inequality holds true, the stationary solution u* is linearly asymptot- 
ically stable. 



6. Direct approach: dissipativity 

Since in the general case of environment dependent vital rates the charac- 
teristic equation is not explicitly available, we shall pursue a different path to 
obtain asymptotic stability. This approach is based on dissipativity calculations 
in the underlying state space X = L 1 (0,m) and proceeds parallel to similar de- 
velopments for the linear semigroup of fiber spinning in Tfl 118) . An added 
advantage of this technique is that we can discuss linear stability of the trivial 
stationary solution and that we can forego imposing positivity conditions on 
the semigroup. In addition we can include the case a = 1 without technical 
difficulties. To our knowledge dissipativity estimates have so far not been used 
in the case of hierarchical size-structured population models. 
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Theorem 6.1 A stationary solution u* is linearly asymptotically stable if 
(*(s,Q*(s)) > w(s) ((1 - a)7 Q (s,Q*(s))tt*(s) + IMQ) 



+ 7(0,Q.(0)) 



P Q {r,Q*{r))u*{r)dr (6.1) 



+ w(s) / P Q (r, Q^{r))u^{r)dr 
Jo 



for < s < m. 

As before the norm on L 1 (0,m) is denoted by || • ||. 

Proof. We will show that, under the given condition, there exists n > such 
that the operator A + B + C + kI is dissipative. Consequently, the semigroup 
{T(t)} t >o generated by the operator A + B + C obeys 

||T(i)|| <e- Kt , t>0, (6.2) 

which proves the claim. 

To obtain dissipativity, assume that, for given h E X, v E Dom(*4) is such 
that, for some A > 0, 

v-X(A + B + C + Kl)v = h. (6.3) 

Then we have 

f-m 

||t>|| = / v(s) sgnv(s) ds 
Jo 

f-m f-m 

= I h{s) sgnu(s) ds - A / 7(3, Q*(s)) v s (s) sgnw(s) ds 
Jo Jo 

f-m 

— A / p*(s) v(s) sgnw(s) ds 
Jo 

-Ay cr*(s) J w(rj) v(r/) di] (6.4) 

/m \ f-m 

w{rj)v{rj) di] \ sgnv(s)ds + \n / v (s) sgnv (s) ds. 

Here we have used the definition sgnO = 0. The set of points in the interval 
(0, m) where v is nonzero is the countable union of disjoint open intervals (di, hi) 
on each of which either v > or v < holds true such that v(di) = for all 
i unless a» = 0, and such that v(bi) = unless bi — in. If (aj,6j) is any such 
interval on which v > we have after integration by parts 

rbi 

v(s)ds< / \h(s)\ds- A7(6 i ,Q*(6 i ))u(6 i ) + A7(a i ,Q»(a i ))t)(o i ) 
J di 

+ A / (k - fi(s, Q*(s)) - (a - 1) w(s)7q(s, Q*(s))u*(s)) u(s)ds (6.5) 
+ A / |cr*(s)|rfs / w(s) |u(s)| rfs. 
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Similarly, on any interval (eii, bi) where v < we have 

bi rbi 

\v(s)\ds< / \h(s)\ds + \^(bi,Q*(bi))v(bi) — Xj(a,i,Q*(a,i))v(ai) 



+ A / (k - /x(s,Q*(s)) - (a - 1) iy(s) 7q(s,Q*(s)) m*(s)) |w(s)| ds (6.6) 



+ A / |o-*(s)| ds / u;(s) |u(s)| ds. 

Finally, noting that u(ai) = = i>(frj) unless = 0, 6j = m, we combine these 
two estimates to obtain 



IMI < |H| + a 7 (o,q,(o))Ko)| 



(6.7) 



+ A / (/t - (j,(s,Q*(s)) + w(s) ((1 - a)7g(s, Q*(s))u*(s) + ||tr*||)) |«(s)|ds. 
Jo 

Since 

|«(0)| = |A(«)| < 



/3(s,Q*(s)) + aw(s) / 0Q(r,Q*(r))u*(r)dr 
o 



+(1 - a) iu(s) / 0Q(r,Q*(r))u*(r)dr 
Jo 



\v{s)\ds (6.8) 



and since condition (|6.1[) is satisfied, we can choose k > such that, for < 

s < m, 

k - fi(s, Q*(s)) + w(s) ((1 - a) 7q(s, Q*(s)) u*(s) + ||ct*||)) 



7(0,Q*(0)) 



+ (1 - a) to(s) / Q {r,Q*(r))u*(r)dr 



P(s,Q*(s)) + aw(s) l3 Q (r,Q4r))u*(r)dr (6.9) 



< 0. 



For such k, we have the desired inequality 

\\v\\<\\h\\ forA>0, 



(6.10) 



thus establishing dissipativity. Hence we conclude that A+B + C + kI generates 
a contraction semigroup. 

□ 



Remark 6.2 Suppose u* is a stationary solution such that condition (|6.1[) is 
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satisfied with (3q > 0. Then we have 

/3(s,Q*( S )) 7 (0,Q*(0)) 



R{Q*)= / 

Jo 



< 



(i{s,Q*(s)) 



exp 



/i(r,Q*(r)) 



o 7( s ,Q*(s)) ° xp V l( r ,Q*( r )) 



/i(r,Q*(r)) 



= 1 — exp — 



dr j ds 
di i ds (6.11) 

r/r | i . (6.12) 



Jo l{r,Q*(r)) 
Hence in light of (|2.2[) we have to conclude that li* = 0. 

Remark 6.3 For the stability of the trivial equilibrium = the criterion 
(16. ip reduces to 

/i(s,Q) > 7(0,0) £(s,0), se[0,m]. (6.13) 

Note that (|6. 13[) clearly implies P(0) < 1, which is the well-known stability 
criterion of the trivial steady state in scramble competition, see |19j . 

Remark 6.4 In scramble competition (a — 1) the stability criterion (|6.ip for 
a stationary solution u* with total (weighted) population 



P* = / w{ri) u*(r)) dr/ 



reads 

/i(s, P*) > w(s) 11(7*11 
where 



0(s,P*)+w(s) / ^ P (r,P»)u*(r)dr 



£(s,P)= 7 (0,P)/3(s,P). 



(6.14) 

< s < m, 
(6.15) 

(6.16) 



Example 6.5 We will give a nontrivial example of a stationary solution for 
which the stability criterion (|6.ip holds true. We choose 



m = l, a = \i w=l, 7(s) = l-is, /i=l 
and let (3 e C ll ([0, 1] x [0, oo)) be positive such that 

Q) = ?S C 1 + s ) ( 2 - 2 Q) if Q < ?, < a < 1. 
159 4 

Then we have again the stationary solution tt*(s) = 1 — o s with 

* 2 « 3 3 



Q*(s) 



< - for < s < 1. 
2 4 ~ 4 ~ ~ 



(6.17) 



(6.18) 



(6.19) 
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Now, however, the positivity condition (|3.13[) is violated. Nonetheless we obtain 



0{s,Q*(s)) + - I Q (r,Q*(r))u*(r)dr+ / Q (r, Q,(r)) u,(r) dr 



160 
159 



s 3 s 2 s 7 



12 



12 



< 1, 



< s < 1. 



(6.20) 



Hence the stationary solution is linearly asymptotically stable by Theorem 16.11 
A straightforward perturbation argument can be used to extend this example to 
a more complicated situation with environment dependent mortality and growth 
rate. 



7. Conclusion 

In this work we have analyzed the linear asymptotic stability of equilibrium 
solutions of a nonlinear hierarchical size-structured population model. We have 
extended our previous mathematical approach in [T31 HH [H] for the case of 
scramble competition models to the hierarchical case. As conjectured in [TU] for 
general physiologically structured models, we note that the linear asymptotic 
stability of stationary solutions is determined by zeros of a characteristic func- 
tion. When the linear dynamical behavior is governed by a positive semigroup, 
the characteristic function has a dominant real root, unless the spectrum of 
the semigroup generator is empty. As we have seen in Section 4, however, this 
function is not explicitly available (except in special cases). Nevertheless, we 
managed to characterize the spectrum of the linearized operator implicitly, by 
deducing an eigenvalue problem for a second order differential operator. This 
characterization allows in principle to further investigate stability questions by 
numerical techniques in case of concrete model ingredients. 

To overcome the severe limitations caused by the spectral characterizations 
of asymptotic stability, we have given a direct dissipativity condition in terms of 
the model ingredients in the relevant state space, guaranteeing the exponential 
decay of the governing linear semigroup. This elementary, though important cri- 
terion allows us to expand linear stability studies beyond the setting of positive 
semigroups and dominant eigenvalues. 
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